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ABSTRACT 


'A  plane-stress  analysis  of  orthotropic  or  isotropic  beams  is  presented. 
The  loading  conditions  considered  are:  (1)  a  concentrated  normal  load 
arbitrarily  located  on  the  beam,  and  (2)  a  distributed  normal  load 
covering  an  arbitrary  length  of  the  beam.  Other  loading  conditions  can 
similarly  be  considered.  Numerical  results  based  on  the  present  analysis 
exhibit  close  agreement  with  existing  experimental  data  from  Sitka 
spruce  beams. 
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INTRODUCTION 


The  solution  of  the  plane-stress  or  plane-strain  problem  for  an  infinite 

3/ 

strip  of  orthotropic  material  was  studied  by  Green  (7).-  In  his  solu¬ 
tion,  both  the  upper  and  lower  boundaries  of  the  strip  must  be  loaded,  a 
situation  not  commonly  encountered  in  engineering  applications.  Using  a 
similar  approach  by  Lamb  (10)  for  a  narrow  beam  of  isotropic  material, 
Smith  and  Voss  (12)  solved  the  problem  of  a  beam  of  either  isotropic  or 
orthotropic  material  with  a  concentrated  load  applied  at  the  center  of 
the  span.  Conway  (2)  used  polynomial  stress  functions  in  his  solutions 
of  an  orthotropic  cantilever  beam  loaded  at  the  free  end  and  a  simply 
supported  beam  under  uniformly  distributed  loads.  He  also  solved  the 
problem  of  a  deep  beam  of  orthotropic  material  symmetrically  loaded  with 
respect  to  the  beam  axis  (3).  Functions  were  derived  by  Silverman  (11) 
to  describe  the  stresses  and  displacements  of  orthotropic  beams  under 
distributed  polynomial  loads.  Hashin  (8)  developed  a  method  to  solve 
the  elasticity  problem  of  long  plane  anisotropic  rectangles  with  continu¬ 
ous  polynomial  stresses  prescribed  on  the  long  sides,  and  force  and 
moment  resultants  prescribed  at  the  ends.  The  problems  analyzed  in 
(2,11)  were  discussed  in  detail  for  anisotropic  material  in  (8).  The 
finite  element  technique  was  applied  by  Hooley  and  Hibbert  (9)  to  solve 
the  plane-stress  problem  of  a  wood  beam.  This  numerical  technique, 
while  applicable  to  intricate  boundary  conditions,  is  limited  by  the 
computer  capacity;  also,  in  areas  of  high  stress  concentrations  where 
interpolation  procedures  must  be  followed  for  stress  estimations,  the 
method  cannot  yield  accurate  results. 
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For  narrow,  rectangular  beams  of  isotropic  material,  solutions  for 
some  specified  loading  and  support  conditions  can  be  found  in  (13). 
Donnell  (5)  applied  infinite  series  for  stresses  in  his  solution  with 
the  first  terms  corrresponding  to  elementary  theory  and  the  latter  terms 
to  increasingly  minor  refinements.  However,  his  solution  cannot  be 
applied  where  discontinuities  exist,  and  the  most  serious  such  discon¬ 
tinuity  is  that  caused  by  a  concentrated  load.  Boley  and  Tolins  (1) 
used  an  iterative  procedure  to  analyze  the  two-dimensional  beam  problem. 
The  applied  loads  considered  by  them  may  consist  of  either  normal  or 
shear  forces  varying  smoothly  along  the  span. 

Experiments  were  performed  by  Cowan  (4)  to  determine  the  horizontal 
shear  stress  distribution  in  beams  of  Sitka  spruce.  He  found  that  in 
the  vicinity  of  a  bearing  support,  the  shear  stress  distributions  were 
considerably  different  from  what  the  simple  beam  theory  would  predict. 
The  same  phenomenon  was  observed  in  (12)  to  exist  in  the  vicinity  of  a 
load  point. 

In  the  study  presented  herein  the  practical  loading  conditions  not 
discussed  in  the  cited  literature  will  be  considered,  namely:  (1)  a 
concentrated  normal  load  arbitrarily  located  on  the  beam,  and  (2)  a 
distributed  normal  load  covering  an  arbitrary  length  of  the  beam.  With 
these  solutions  available,  a  large  number  of  practical  beam  problems  can 
be  handled  using  the  method  of  superposition. 

THEORETICAL  ANALYSIS 


The  beam  considered  is  assumed  to  be  an  orthotropic  solid  having 
its  edges  parallel  to  two  perpendicular  axes  of  elastic  symmetry.  The 
thickness  of  the  beam  is  assumed  small  as  compared  with  the  vertical 
depth  so  that  the  problem  can  be  treated  as  one  of  plane  stress. 

The  x-axis  is  taken  along  the  center  line  of  the  beam  and  the 

y-axis  coincident  with  the  left  end  of  the  beam  as  shown  in  figure  1(a). 

For  the  isotropic  case,  the  pertinent  derivations  are  presented  in 
Appendix  I. 

STRESS  AND  LOAD  FUNCTIONS 


For  the  state  of  plane  stress  in  the  orthotropic  beam,  the  stress 
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Equation  1  may  be  satisfied  by  taking  the  function  <(>  in  the  form 


<f»  =  ?  f(n)  sin  (2) 

n=l  x 


In  which  f(q)  is  a  function  of  q  only.  Substituting  equation  2  into 
equation  1  and  using  the  notation  m  =  n n/i,  one  finds  the  following 
equation  for  determining  f(q): 


m‘f  -  2K.2  ^  ^  ,  0 

dn2  dq4 


(3) 


The  general  integral  of  this  linear  differential  equation  with  constant 
coefficients  is 


f(H)  =  Ancosh  maq  +  Bncosh  mpq 

(4) 

+  Cnsinh  maq  +  Dnsinh  mpq 
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where 


and 


(4a) 


(4b) 


The  stress  function  is  then  obtained  by  substituting  equation  4  into 
equation  2,  and  the  corresponding  stress  components  are 
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in  which  the  constants  A  ,  B  ,  C  ,  and  D  are  to  be  determined  from  the 

n’  n*  n*  n 

loading  conditions  and  t  represents  x  . 

Assuming  that  the  load  function  can  be  expressed  as  a  Fourier  series, 


p(x)  =  ?  ansin  nr  (g) 

n=l 


the  coefficients  an  in  equation  8  can  be  written  as 


a 

n 


(8a) 


Let  p  (x)  be  represented  by  P/2 £  acting  at  x  =  a  as  shown  in  figure  1(b) 
Then 


a 

n 


(8b) 


For  a  concentrated  load  P  acting  at  x  =  a,  one  obtains  from  equation  8b 


i  ■  2P  nna  . 

lim  a  =  y-  sin  (8c) 

1+ 0  £  £ 


Hence,  for  any  specified  loading  condition,  the  load  function  can  be 
expressed  as  a  Fourier  series  following  the  above  procedure  to  be  consis 
tent  with  the  stress  function  defined  by  equation  2. 
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METHOD  OF  SOLUTION 


To  determine  Aq,  Bq,  Cq,  and  Dq  in  equations  5,  6,  and  7, 

four  boundary  conditions  are  needed.  Consider  the  problem  described  in 
figure  1(b).  The  four  boundary  conditions  are: 


(1) 

Oy  =  P(x) 

at  y  =  h 

(9) 

(2) 

o 

II 

at  y  =  h 

(10) 

(3) 

a  =  0 
y 

at  y  =  -h 

(ID 

W 

T  =  0 

at  y  =  -h 

(12) 

From  these  boundary  conditions,  one  obtains 


(13) 


where 

D  =  m2[l  -  ^(a*-  +  P^)]  cosh  2meh(a+p) 
+m2[l  +  ^(a2  +  02)]  cosh  2meh(a-0) 


-  2m2  (13a) 

a  , 

[(0  -l)cosh  meh(a+20) 

-  (02+l)cosh  meh(a-20)  +  2cosh  math]  (13b) 


Db  =  (a2-l)cosh  meh(2o+0) 
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-  (a  +l)cosh  meh(2a-0)  +  2cosh  mpeh]  (13c) 


a  , 

Dc  =  (P  -l)sinh  meh(a+2p) 

2 

-  (P  +l)sinh  meh(a-2p)  -  2sinh  roaeh]  (13d) 
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D  =  — [(a-l)sinh  m£h(2a+p) 
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+  (a  +l)sinh  meh(2a-p)  -2sinh  mpeh]  (13e) 


It  is  noted  that  equations  13  are  applicable  for  an  corresponding  to  any 

load  function.  It  is  also  noted  that  if  the  load  were  applied  on  the 
lower  side  of  the  beam,  all  that  need  be  done  is  to  change  the  sign  for 
h  in  equations  9-12. 

The  infinite  series  obtained  by  substituting  equation  13  into 
equations  5,  6,  and  7  can  easily  be  solved  using  a  digital  computer.  In 
order  to  save  on  computer  storage  and  time,  a  method  of  reducing  the 
infinite  series  to  the  sum  of  a  finite  series  and  a  closed  form  is 
described  in  Appendix  II. 


DEFLECTION  OF  BEAM 


Deflection  of  a  beam  may  indicate  its  performance  under  specified 
loading  conditions.  It  may  also  serve  to  derive  the  reactions  in  an 
indeterminate  situation.  From  the  strain-displacement  relationships  and 
the  generalized  Hooke's  law,  the  following  equations  can  be  obtained: 

V  =  |  (ax  -I>,xyay)dx  +  r(y)  (14) 

Eyv  =  /  (ay  -^yx°x)dy  +  s(x)  (15) 


where  u  and  v  are  the  displacements  in  the  x-  and  y-directions , 
respectively.  Also, 
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From  the  stress  expressions  in  equations  5,  6,  and  7 
equation  16  the  following  relation  must  hold  true 
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Likewise, 
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s  =  -E  K  x  +  K„  (19) 

y  o  2 

The  constants  Kq,  K^,  and  K2  in  equations  18  and  19  are  to  be  determined 
from  the  known  displacements  at  the  supports. 

Usually,  only  the  vertical  deflection  of  the  beam  center  line  is  of 
interest  in  engineering  applications.  One  then  obtains  from  equations  5, 
6,  15,  and  19,  the  following  expression 


E  v  =  -t  m[C  (—  + 1/  ea)  +  D 
y  j  1  n  ea  yx  nv£0 


+V  eB)]*sin  mx  -  E  K  x  +  K„ 
yx  v  y  o  2 


(20) 


For  a  beam  whose  vertical  displacements  at  the  two  ends  are  zero,  it  is 
seen  that  both  Kq  and  in  equation  20  must  vanish. 


NUMERICAL  RESULTS 


Numerical  calculations  were  made  for  the  wood  beams  used  in  Cowan's 
experiments  (4)  for  shear  stress  measurements  using  a  Univac  1100  computer. 
In  his  experiments,  a  simple  beam  and  a  two-span  continuous  beam  were 
used,  which  were  clear  straight-grained  specimens  of  Sitka  spruce. 

Supports  and  load  blocks  also  were  of  Sitka  spruce  to  simulate  actual 
construction  conditions. 

For  the  simple  beam,  the  applied  load  and  reactions  and  the  geo¬ 
metrical  dimensions  from  (4)  are  shown  in  figure  2.  The  right  support 
is  a  roller-type  support,  but  the  left  support  is  a  flat  bearing  block 
restrained  against  horizontal  movement.  Hence,  the  left  support  could 
be  assumed  to  approximate  a  hinged  support  with  resultant  of  the  reaction 
forces  acting  closer  to  the  load  point.  The  reaction  forces  at  the  left 
support  are  therefore  assumed  to  be  distributed  as  shown  in  figure  2(b). 

The  material  properties  taken  from  (4)  were  as  follows: 


9 


1.84  x  106  psi  (1.2686  x  1010  N/m2) 
6.62  x  10^  psi  (4.5643  x  10®  N/m2) 
1.089  x  105  psi  (7.5083  x  108  N/m2) 
0.26 
0.0094 


the  orthogonal  shear  stress  distributions  in  the  beam 
at  the  left  support.  The  theoretical  predictions  and  the  experimental 
data  are  seen  to  agree  closely  with  each  other.  The  differences  are, 
of  course,  due  to  the  material  heterogeneity  and  the  assumed  stress 
distribution  at  the  support.  It  is  to  be  noted  that  the  stresses  in 
the  beam  are  strongly  dependent  on  the  stress  distribution  at  the  support 
which  can  only  be  estimated. 

The  vertical  compressive  stress  distributions  are  shown  in  figure  4. 

By  comparing  figures  3(c)  and  3(d)  with  figure  4,  it  is  seen  that  the 
large  shear  stresses  do  not  always  occur  where  the  vertical  compressive 
stresses  are  high. 

The  bending  stresses  plotted  in  figure  5  are  seen  to  be  highly  non¬ 
linear  and  are  negative  in  both  the  upper  and  lower  regions  of  the  beam 
above  the  support. 

Figure  6  shows  the  applied  loads  and  reactions  and  the  dimensions  of 
the  two-span  continuous  beam.  The  middle  reaction  is  a  flat  bearing 
block  of  wood  8  inches  long.  However,  due  to  material  variabilities, 
the  resultant  of  the  reaction  forces  was  found  to  be  2  inches  on  the  right 
of  the  center  of  the  middle  support.  In  the  numerical  calculations,  the 
reaction  forces  are  then  assumed  to  be  distributed  as  described  in 
figure  6(b) . 

The  material  properties  for  this  beam  were  the  same  as  those  for  the 
simple  beam  except  the  modulus  of  rigidity  was  found  to  be 

1.121  x  105  psi  (7.7289  x  108  N/m2)  rather  than  1.089  x  105  psi 

(7.5083  x  108  N/m2)  according  to  (4). 

The  shear  stress  distributions  at  the  middle  support  are  shown  in 
figure  7.  Due  to  the  material  heterogeneity  and  the  inherent  experimental 
inexactness,  considerable  randomness  exists  in  the  experimental  data. 
However,  the  variations  of  the  experimental  shear  stress  with  location 
follow  generally  the  same  trend  as  predicted  theoretically. 
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The  vertical  compressive  stress  distributions  plotted  in  figure  8 
disclose  their  strong  dependency  on  the  assumed  reactive  forces  at  the 
support.  This  is,  of  course,  partly  due  to  the  relatively  large  distances 
of  the  applied  loads  and  the  end  supports  from  the  middle  support. 

At  a  given  depth  of  the  beam,  the  nonlinear  bending  stresses  are 

seen  to  vary  only  slightly  across  the  length  of  the  support  as  described 

in  figure  9. 

Numerical  calculations  were  also  performed  for  the  simple  beam  of 
Sitka  spruce  loaded  at  the  center  as  reported  in  (12) .  For  a  concentrated 
central  load,  the  theoretical  predictions  of  the  shear  stresses  in  the 
vicinity  of  the  load  point  of  the  present  study  and  of  reference  (12) 
are  essentially  the  same.  In  (12),  the  following  function  was  used  for 
a  unit  load  acting  on  the  beam: 


0(x) 


7t  2.  e2 
x  +  o 


(21) 


where  x  is  measured  from  the  center  of  an  infinitely  long  beam  along  its 
axis  and  the  value  of  6  is  related  to  the  length  covered  by  the  load. 

For  6=0,  the  load  is  concentrated;  for  small  values  of  6,  equation  21 
was  assumed  suitable  for  representation  of  a  load  applied  over  a  small 
area  by  a  curved  loading  block.  In  the  test  reported  in  (12),  the  central 
load  was  applied  over  a  length  of  about  0.4  inch  (10.16  mm)  by  a  cylindrical 
block.  At  a  point  close  to  the  load,  the  calculated  shear  for  6=0  inch 

(0  mm)  wa*  235  psi  (1.6202  x  10^  N/m^)  and  for  6  =  0.25  inch  (6.35  mm) 

5  2 

was  116  psi  (7.9978  x  10  N/m  ),  which  was  close  to  the  measured  value. 

Based  on  the  present  work,  the  calculated  shear  for  £  =  0  inch  (0  mm)  was 

234  psi  (1.6134  x  10^  N/m^),  for  £  =  0.2  inch  (5.08  mm)  was  229  psi 

(1.5789  x  10^  N/m^)  and  for  £  =  0-3  inch  (7.62  mm)  was  212  psi 
6  2 

(1.4617  x  10  N/m  ).  To  reduce  the  shear  stress  to  the  measured  value, 
the  length  covered  by  the  load  must  be  much  larger  than  0.6  inch 
(15.24  mm).  It  is  believed  that  the  measured  value  at  that  point  is 

likely  to  be  in  error  and  6  cannot  be  properly  used  to  represent  a  load 

over  a  specified  small  length.  In  fact,  the  statement  made  in  (12) 
that  84  percent  of  the  area  under  the  curve  of  equation  21  is  over  a 
length  26  on  either  side  of  x  =  0  is  wrong.  It  can  be  shown  that  for  a 

length  of  26,  the  percentage  is  50;  for  a  length  of  46,  the  percentage 

is  70;  and  for  a  length  of  86,  the  percentage  is  84. 
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For  the  isotropic  case,  the  calculated  stresses  and  deflections 
based  on  the  present  work  agree  with  the  best  estimations  in  (13)  for  a 
simple  beam  loaded  at  the  middle. 


DISCUSSION  AND  CONCLUSION 


The  present  analysis  demonstrates  a  procedure  for  calculating  the 
stresses  and  deflection  of  a  two-dimensional  orthotropic  or  isotropic 
beam  under  arbitrary  loading  conditions.  In  the  solutions,  each  of  the 
infinite  series  has  been  reduced  to  the  sum  of  a  finite  series  and  a 
closed  form,  making  it  possible  to  obtain  accurate  numerical  results 
using  a  computer  of  relatively  small  capacity. 

Numerical  results  based  on  the  present  work  indicate  that  they  are 
in  reasonable  agreement  with  existing  experimental  data.  For  specified 
loading  conditions  that  can  be  considered  by  existing  theories,  e.g., 
(12)  and  (13),  the  present  method  can  yield  the  same  predictions  as  the 
other  theories. 

The  method  is  very  useful  in  predicting  the  stress  distributions  in 
a  beam,  especially  in  the  vicinity  of  a  support  or  a  nominal  load  point. 
However,  for  the  latter  purpose,  the  reactive  stresses  at  a  support  or 
the  acting  stresses  at  a  nominal  load  point  must  be  correctly  estimated 
in  order  to  obtain  accurate  results. 
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Figure  l.--(a)  Coordinate  system  and  geometry  of  beam 
(b)  load  function  and  load  location  on  beam. 
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790.8  LB 


766.8  LB 


4634" 

48" 

(a) 


Figure  2. --Simple  beam  (a)  applied  load  and  reactions;  (b) 
assumed  compressive  forces  at  left  support  (1  in.  =  25.4  mm; 
1  lb.  =  4.4482  N;  1  lb. /in.  =  175  N/m). 


(M  146  767) 


(M  146  768,  146  769,  146  770,  146  771) 


DEPTH,  IN. 


(c) 


Figure  3. --Variations  of  shear  stress  vs.  depth 
from  top  of  simple  beam  (1  in.  =25.4  mm; 

1  psi  =  6.8947  x  103  N/m2).  (Continued) 

(M  146  768,  146  769,  146  770,  146  771) 
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arbitrarily  located  on  the  beam,  and  (2)  a  distributed  arbitrarily  located  on  the  beam,  and  (2)  a  distributed 

normal  load  covering  an  arbitrary  length  of  the  beam.  normal  load  covering  an  arbitrary  length  of  the  beam. 


6  7  8  9  10  II  12 

DISTANCE,  IN. 

Figure  4. — Theoretical  predictions  of  vertical  compressive 
stress  vs.  distance  from  left  end  of  simple  beam;  at  left 
end  of  support  x  =  4.46  inches;  at  right  end  of  support 
x  =  10.46  inches.  (1  in.  =  25.4  mm;  1  psi  =  6.8947  x  103  N/m2). 

(M  146  772) 
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Figure  5 . --Theoretical  predictions  of  bending  stress 
vs.  depth  from  top  of  simple  beam  (1  in.  =  25.4  mm; 
1  psi  =  6.8947  x  103  N/m2!. 
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(1  in.  =  25.4  mm;  1  lb.  =  4.4482  N;  1  lb. /in.  =  175  N/m) . 
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Figure  7 . --Variations  of  shear  stress  vs.  distance 
from  left  end  of  continuous  beam  (1  in.  =  25.4  mm; 
1  psi  =  6.894/  x  103  N/m2).  (Continued) 

(M  146  775,  1.46  776,  146  777,  146  778) 


23 


s  tress,  psi 


DISTANCE  FROM 
TOP  OF  BEAM,  IN 

SHEAR  STI 
EXPERIMENTAL 

’ESS,  PS! 
THEORETICAL 

5 

o 

6 

A 

— 

80 


Figure  7 . -“Variations  of  shear  stress  vs.  distance 


from  left  end  of  continuous  beam  (1  in.  =  25.4  mm; 
1  psi  =  6.8947  x  103  N/m2).  Continued) 

(M  146  775,  146  776,  146  777,  146  778) 
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DISTANCE  FROM  TOP 


Figure  7 .--Variations  of  shear  stress  vs.  distance 
from  left  end  of  continuous  beam  (1  in.  =  25.4  mm; 
1  psi  =  6.8947  x  103  N/m2).  (Continued) 

(M  146  775,  146  776,  146  777,  146  778) 
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5  TRESS,  PS  I 


Figure  8. — Theoretical  predictions  of  vertical  compressive 
stress  vs.  distance  from  left  end  of  continuous  beam;  at 
left  end  of  middle  support  x  =  77.5  in.;  at  right  end  of 
middle  support  x  =  85.5  in.  (1  in.  =  25. A  mm;  1  psi  = 
6.8947  x  10*  N/m2). 

(M  146  779) 
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Figure  9 . --Theoretical  predictions  of  bending  stress 

vs.  depth  from  top  of  continuous  beam  (1  in.  =  25.4  mm 
1  psi  =  6.8947  x  103  N/m2). 
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APPENDIX  I . —FORMULATION  OF  ISOTROPIC  BEAM 


For  an  isotropic  beam,  one  cannot  obtain  a  solution  by  simply  setting 

a=P=£=lin  equations  13  and  then  substitute  them  into  equations  5, 

6,  and  7,  since  this  will  cause  the  constants  A  ,  B  ,  C  ,  and  D_  in  these 
’  ’  n  n  n  n 

equations  to  become  indeterminate.  To  effect  the  conversion,  one  must 

set  e  =  1,  P  =  1  and  then  differentiate  D,  Da>  D^,  Dc,  and  in  equations  13 

a 

with  respect  to  a  twice  before  setting  a  -  1.  After  some  tedious  algebraic 
manipulations,  the  following  results  are  obtained: 


o  =  ?  ^ — [D.cosh  my  +  D~(2cosh  my 
X  n=l  D" 


where 


+  mysinh  my)  +  D^sinh  my 
+  D^(2sinh  my  +  mycosh  my)] sin  mx 


(1-1) 


o  =  -?  - — [D  cosh  my  +  D.mysinh  my 
y  n=l  D"  1 

+  D^sinh  my  +  D^mycosh  my] sin  mx 


CI-2) 


T  =  -?  - — [D..sinh  my  +  D.(sinb  my 

n=l  D"  1 

+  mycosh  my)  +  D^cosh  my 


+  D^(cosh  my  +  mysinh  my)]cos  mx 


(1-3) 


D"  =  4m^(8m^h^  -  cosh  Amh  +1) 


(1-4) 
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D.  =  4a  [cosh  3mh  +  mhsinh  3mh 
1  n 

-  (4m^h^+l)cosh  mh  -  3mhsinh  mh]  ( I -5 ) 

=  4an(4mhsinh  mh  +  cosh  mh 

-  cosh  3mh]  (1-6) 

D„  =  4a  [mhcosh  3mh  +  sinh  3mh 
3  n 

2  2 

+  (4m  h  +l)sinh  mh  ♦  3mhcosh  mh]  (1-7) 

=  -4an(4mhcosh  mh  +  sinh  mh  ♦  sinh  3mh]  (1-8) 


in  which  a^  is  expressed  by  equation  8c. 

It  is  to  be  noted  that  if  the  stress  function  for  the  isotropic 
case  as  described  in  reference  (13)  were  used,  one  would  have  obtained 
the  same  results  for  the  same  loading  condition  as  shown  in  figure  lb. 
This  also  serves  to  confirm  that  the  derivations  for  the  orthotropic 
beam  in  the  present  study  are  correct. 

The  vertical  deflection  of  the  center  line  of  the  isotropic  beam 
can  be  shown  to  take  the  following  form: 


mD,  mD, 

Ev  =  2.  [-(1+V) -  +  (147) - ]*sin  mx 

n=l  D"  D" 


-  EK  x  +  K, 

O  2 


(1-9) 


which  corresponds  to  equation  20  for  the  orthotropic  beam. 
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APPENDIX  II.— SIMPLIFIED  FORMS  OF  STRESS  COMPONENTS 


It  is  possible  to  reduce  the  infinite  series  for  each  stress  component 
to  the  sum  of  a  finite  series  and  a  closed  form.  To  do  this,  one  can 
first  express  the  hyperbolic  functions  in  terms  of  exponential  functions 
and  then,  after  some  simplifications,  reduce  the  infinite  series  comprising 
exponential  functions  and  trigonometrical  functions  to  closed  form. 


A.  ORTHOTROPIC  BEAM 


It  can  be  shown  that  equation  13  can  be  reduced  to 


_  an  p2-  1  -meha 
n  m2  2-(a2+p2)6 


(II-l) 


B 

n 


an  a2-  1  -mehp 
m2  2-(a2+p2)  6 


( I 1-2) 


C 

n 


(II-3) 


D  =  B  (II-4) 

n  n 

when  e  2me^  is  negligibly  small  compared  with  unity.  Equations  5,  6, 
and  7  can  now  be  expressed  in  the  following  form  when  a  in  equation  8c 
is  used: 


a 

x 


2 

e  P 

£[2-(ct2+p2)] 


•  t 

n=l 


[a2(p2-l) . 


e-ma(eh-n)  +  p2((|2_1)  .e-mp(£h-q)] . 


[cos  m(a-x)  -  cos  m(a+x)) 


(II _5) 


( 1 1  “6 ) 


O  = - ?  K^-De-^-n) 

y  SL l2-(ci  +pz)]  n=l 

+  (a2-l)e'mp(eh-n)]-(cos  m(a-x)- 
cos  m(a+x)J 

- - ?  l«02-l)e'IMl£h-,» 

£|2-(a *p‘)l  n=l 

+  0(a2-l)e  r^]»[sin  m(a-x) 

+  sin  m(a+x)]  ( I I - 7 ) 

For  other  expressions  of  a  ,  the  above  equations  need  to  be  modified  but 
the  procedures  of  derivation  remain  the  same. 

From  reference  (6),  the  following  identities  can  be  derived: 

f  e”nq6cos  n6  =  Tj+T^  (II-8) 

n=l 

?  e  nq0sin  n6  =  -qT,  +  —  ( I 1-9 ) 

n=l  1  q 

where 


t  -  1  teq  (cos  8  -  qsin  6)  -  1. 

*1  9  l  -.ft  -ft  J  >*1^0 

(l+q  )  eq  (  eq  -  2cos  0)  +  1 


qe  4  fqcos  6  ♦  sin  9  •  q  ,  „ 

2  (l+q2)  e-q0(e-q6-2cose)+l  " 


It  can  be  seen  that  the  series  in  equations  1 1 -5 ,  1 1 -6 ,  and  1 1 - 7  take  the 
same  form  as  either  equation  1 1 —8  or  equation  1 1 -9 -  Therefore,  the  stress 
components  can  ail  be  reduced  to  closed  form,  which,  though  complicated, 
can  easily  be  programmed  for  a  computer. 

It  is  to  be  noted  that  the  first  terms  in  equations  5,  6,  and  7  must  be 

retained  until  n  is  so  large  that  e  2me^  becomes  negligibly  small  when 
compared  with  unity,  and  that  the  same  number  of  terms  must  be  deducted 
from  equations  1 1-5 ,  II-6,  and  1 1 - 7  before  the  closed  form  expressions  are 
added  to  the  finite  series.  All  these  can  easily  be  handled  in  a  computer 
program. 


B.  ISOTROPIC  BEAM 


For  the  isotropic  beam  discussed  in  Appendix  I,  the  following 
results  are  obtained: 


m2D 

- -  =  -a  (mh  +  l)e'mn  (11-10) 

nn  ft 


-mh 

a  e 
n 


(ii-n) 


D3  =  D1 

(H-12) 

=  °2 

(11-13) 

when  e  is  small  compared  with  unity.  Again,  with  a^ 

equation  8c,  the  stress  components  can  be  written  as 

represented  by 

°x  =  nil  l-"(h-y)-l)e"m(h"y) 

•  [cos  m(a-x)  -  cos  m(a+x)] 

(11-14) 
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Oy  =  |  n|j  [m(h-y)+l]e 


(11-15) 


-m(h-y) 

•  [cos  m(a-x)  -  cos  m(a+x)] 

X  =  £  t  m(h-y)e-m(h-y) 

Z  n=l 

•[sin  m(a-x)  +  sin  m(a+x)]  (11-16) 

To  reduce  the  infinite  series  in  these  equations  to  closed  form,  in 
addition  to  equations  II-8  and  1 1 -9 ,  the  following  identities,  also 
derived  from  (6),  will  be  needed: 

t  ne"nq0cos  n0  =  (l-q2)T3  +  2qT^  (11-17) 

n=l 

|  ne'nq0sin  n0  =  -2qT3  +  (l-q2)T4  (11-18) 

where 

T3  =  e  q0  {  (e  2q0-l) [q2(cos  0-1)  +  cos  0]  +  2e  q0(qsin  0-1) 

+  2  (cos  0-qsin  0)}  /  [(l+q2)(e  2q0-2e  q0cos  0  +  l)]2,  q>0 

T  =  eq0  {  (e2q0-l) [ (l-q2)sin  0  +  2qcos  0]  +  4q(cos  0 
4 

-  eq0)}  / [ (1+q2) (e2q0  -  2eq0cos  0  +  l)]2,  q>0 

The  numerical  procedures  are  the  same  as  in  the  orthotropic  case. 
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APPENDIX  I II. --NOTATION 


The  following  symbols  are  used  in  this  paper: 


f(n),  r(y),  s(x) 
G 

xy 

h 


constants  of  integration; 

x-coordinate  of  a  load  point; 

Fourier  coefficients; 

moduli  of  elasticity  in  x-  and  y-directions , 

respectively; 
functions ; 

modulus  of  rigidity  associated  with  xy-plane; 
half  depth  of  beam; 


K  = 


m 

n 

P 

p(x) 
u,  v 

x,  y 


constants ; 

length  of  beam; 
nn/£; 
integer ; 
applied  load; 
load  function; 

displacements  in  x-and  y-directions,  respectively; 
Cartesian  coordinates 
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<t>  =  stress  function; 


Vxy 


=  Poisson's  ratio  associated  with 


stress  in  x-direction  and  strains  in  x- 
and  y-directions ; 


o  a  =  normal  stresses  in  x-  and  y-direc 
x,  y 

tions,  respectively; 


T  =  orthogonal  shear  stress. 
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